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Abstract 

We study the electronic excitations near the charge-transfer gap in insu- 
lating Cu02 planes, starting from a six-band model which includes p-,^ and 
dxy orbitals and Cu-0 nearest-neighbor repulsion Upd- While the low lying 
electronic excitations in the doped system are well described by a modified 
t — J model, the excitonic states of the insulator include hybrid d^y— Pn states 
of A2g symmetry. We also obtain excitons of symmetries Big and E^, and 
eventually Aig, which can be explained within a one-band model. The results 

agree with observed optical absorption and Raman excitations. 
PACS Numbers: 74.25.Gz, 74.72.-h, 71.35.+Z 
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There is general consensus that the main electronic properties of superconducting Cu02 
planes are well described by the three-band Hubbard model ifa^ |I|] or effective one-band 
models derived from it f^-p^. In particular, accurate low-energy reductions of ifsb have 
been made using the cell- perturbation method PJ^JTO[] or fitting energy levels and 
the resulting generalized t — J or Hubbard models were used to explain many fundamental 
properties of these materials, like superconductivity |ll| , angle-resolved photoemission [jl2 



x-ray absorption spectrum in La2-a:Sra;Cu04 fgJT^, metal-insulator transition as a function 
of the charge-transfer energy A P] and Raman scattering |T^ , [T5[| . 



While the general features of optical absorption experiments [16,13] have been explained 



using H^f, [|T8| , recent Raman experiments |T9| in several insulating cuprates show excitations 



below the charge-transfer gap which cannot be explained with H^f,. This raises serious doubts 
on the validity of the very large amount of theoretical work ultimately based on this model. 
The most intense observed Raman transitions have symmetries and Big . The former 
lies lower in energy, is more intense for the frequency of the incident light used, and has 
been ascribed to transitions from dx2_y2 to d^y or orbitals ||TP[. The latter two orbitals 
are not contained in H^f,. The Big peak follows closely the optical absorption edge and its 
origin is not clear so far. 

The purpose of this work is twofold. First, to explain within a unifying approach the 
Eu exciton observed in optical absorption experiments, and the main features of the Raman 
experiments in insulating systems. Second, to analyze the more fundamental issue of how to 
reconcile the A2g peak with our present theoretical understanding of the cuprates. We study 
a six-band model Hqi, obtained generalizing H31, to include the dxy and Pt, orbitals. The model 
is treated within the cell-perturbation method 0]^,^. This method has been shown to give 
accurate results and at the present, we believe that there are no accurate alternative 
methods to handle this particular problem. We obtain in agreement with experiment that 
the observed absorption Eu peak and the Raman Big peak lie very close in energy, and 
(0.1-0.2) eV below the absorption edge. While these peaks could be explained within an 
effective generalized one-band Hubbard model, the Raman A2g peak corresponds to an 



excitation of a hole from a (predominantly) dx2_y2 orbital to a linear combination of d^y 
and holes with symmetry. However, holes added to the ground state of the insulating 
system enter predominantly 2pcr orbitals in agreement with experiment and our previous 
understanding of these systems. 

The six-band Hamiltonian can be written as: 



Heb = ^xy^iaa + XI {^^jp<T + A^rijva) + Ud ^ip^ni pi^' + Up ^ 

icr ja i,l3(Tyt (3' a' j,aa^a'a' 

+ Upd ^il3ani+Saa' + tpd^ f (pUsadia + h.C.) + 9{^){A+5aO'i<T + h.C.) 

iSa/Saa' i§cr iSa 

+ J2H^)(^PpP]+l^Pj^ + " tp^{n]^^^Pj^ + h.C.)] (l) 

with ni0^ = Ptpia, p = doT a; rijaa = oi]„aj^, a = p om; dl„, aj^, p]^,7rj^ create a 3d^2_y2, 
3dxy, 2po- or 2pt, hole respectively with spin a at Cu site i or O site j. The indices i + 6 (j + 7) 
label the four O atoms nearest to the Cu (O) atom at site i (j). The functions f{S), g{6) 
and h{6) describe the sign of the hoppings as functions of direction. 

Most of the parameters of the model are known, at least approximately, from constrained- 
local-density-functional calculations P PT|P^ . The parameters involving the orbitals are, 
according to the study of Mattheis and Hamann [^, tp^ = 0.45tpp, tTrn = 0.48 tpp, and 

= A — leV. We take tpd — l.2eV as the unit of energy and the following parameters 
iJII-ll: A = 2.35, A,^y = 1.35, A^ = 1.55, Ud = 8.33, Up = 3.33, Upd = 0.833, Ua = 



0.65, tpp 0.5, ^7^77 0.3, tpTY t-jj-jj-. 

The large number of orbitals per unit cell and the highly correlated nature of H^f,, 
make impossible at present an exact treatment of the model in a large enough system . 
However, the methods we use are accurate enough to explain the nature of the observed 
optical excitations. Our basic approach is the cell-perturbation method PJ^,^. For each 
Cu site i, four orthogonal Wannier functions of the O 2p orbitals centered at that site are 
defined (see Fig. 1). Two of them, of symmetry big (like the 3(i^.2_y2 orbital) and aig ("non- 
bonding") correspond to the p„ orbitals. The Wannier functions constructed with the p,^ 
orbitals have symmetry a2g and 629- The cell Hamiltonian Hi contains all interactions and 
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hoppings involving only the six orthogonal Wannier functions centered at site i. For each 
cell, Hi is diagonalized exactly and the intercell part is written in terms of the eigenstates 
of lowest energy of Hi. The advantage of this method is that Hi contains the highest energy 
in the problem {Ud) and most of the Upd and hopping terms. For one hole in the cell, the 
ground state of Hi, \igcr) has predominantly Cu character and is a linear combination of 
both big states with spin a The first excited state liner) corresponds to one hole 

in a linear combination of the d^y and O Wannier functions and lies 1.63 above {igcr) 
for the parameters given before. For two holes, the ground state of Hi, \ig2) is a linear 
combination of the three singlets constructed with both big orbitals, and corresponds to the 



Zhang-Rice singlet (ZRS) of the three-band model and one-band models derived from 
it The first excited state lying 0.99 above, denoted as |i7r2), is spin degenerate 

and contains one a2g hole and one big hole of predominantly Cu character. The hopping 
\ig2) \jgl) \w^) |j7r2) vanishes by symmetry if i and j are nearest neighbors (NN). The 
larger hopping of this kind is for next-NN and is OAbtp^^ = 0.13, considerably smaller than 
the energy difference 0.99 and about a third of the hopping between ZRS. Thus, added 
holes in the system form ZRS and the amount of holes in the ground state of hole-doped 
systems is very small for small doping. The same happens for d^y holes. For the explicit 
calculation of the exciton energies we neglect these small hoppings. 

The \ig(j) , \ig2) states are mapped into a generalized one-band model Hn, which includes 
NN, next-NN hopping and NN repulsions which depend on the occupation of the sites 
involved The NN repulsions are of order Upd/ 4:. The movement of a vacant site (Cu^ 

state) and that of a ZRS in the insulating system are described by corresponding t — t' — J 
models with J = 0.12 (near to the experimental value) and the other parameters taken from 
ifif,: for the movement of the ZRS (vacant site) we obtain = 0.39, = 0.034 {t^ = 
0.34, = 0.005) in units of tpd- The problem of the exciton energy is then approximated by 
that of an electron and a hole moving freely in their respective quasiparticle bands, except 
at NN, where they feel an attraction Vi ~ 0.18 and an on-site interaction determined by the 
eigenstates of Hi. A similar approach has been used previously by Belinicher et al. [Q. The 



quasiparticle band is taken of the form p6|,|T2|: 



^ki = Coi + 4^2* COS COS ky + 2t3j(cos 2kx + cos 2ky) (2) 

where i = h{e) for holes (electrons). The parameters are determined from t, t' and J by the 
self-consistent Born approximation. We obtain e^h = —0.507, = 0.0085, t^h = 0.017 
(eoe = —0.507, t2e = 0.021, t-^e = 0.015). The resulting energy gap is 1.39. 

For total wave vector K = 0, the problem takes the form of a one-particle one. For 
symmetry Aig we do not obtain a bound state because Vi is slightly smaller than the critical 
value (0.19) necessary to obtain an exciton. For symmetries Big and Eu the resulting 
binding energies are 0.09 and 0.11 respectively. A scheme of the quasiparticle bands and 
these excitonic energies is represented at the right of Fig. 2. It can be shown that the 
binding energies of the E^, Big and Aig excitons are mainly related with the quasiparticle 
energies at ('/r/2, 7r/2), (vr, 0) and (0, 0) respectively. The lower the quasiparticle energy of 
the t — t' — J model at one of these points, the higher the binding energy of the respective 
mode. 

The amplitude of the excitonic wave functions for each irreducible representation is shown 
in Fig. 3. For symmetry Aig, there is a small probability of finding the electron and hole at 
the same site, in the state \igcr > or an excited state of Hi of symmetry big and energy ~ 4 
above \iga >. Similarly for symmetry Big, there is a small occupation of non-bonding aig 
one-particle states of Hi at energy 4.30 above \iga >. 

The position of the Raman A2g peak (0.18 below the band edge) was calculated in a 
similar way. The main difference is that the holes move in the two bands of states and 
that there is a stronger effective on-site attraction. As a consequence, the binding energy 
of this exciton, with respect to the bands, is larger (see Fig. 2) and the wave function is 
more localized than for the previous solutions ( see Fig. 3). For symmetry B2g, the bound 
state falls above the gap, inside the ZRS band. For our results to be compatible with the 
Raman experiments and the most accepted physical picture of the hole-doped cuprates. 



it is necessary that while the lowest ZRS lies lower in energy than the states containing 



orbitals, the exciton constructed with the latter hes inside the charge-transfer gap. 
This is the case for reasonable parameters, like the ones we have taken. The binding energy 
of the Aig, Big and Eu excitons is only sensitive to Upd and the shape of the quasiparticle 
band. At least in the range 0.4 < Upd < 1.2 both Big and Eu excitons bind and they lie 
less than ~ 0.2 below the gap edge, in agreement with experiments in different compounds 
19| . While the position of the A2g peak is more sensitive to variations in A^yjA^^ or tj^a, 



a different qualitative explanation of it with parameters consistent with those given in the 
literature is unlikely. 

The intensity of the optical conductivity and Raman spectra is determined by matrix 
elements of the kinetic energy in a given direction Tx,Ty, and of the current operator j 



|, |25|JT^ . This operator can be constructed from Hqi, in a standard way p4[ and is linear in 
the hopping terms. An analysis based on the form of the above mentioned operators and the 
eigenstates of Hi, allows to infer the character of the expected observable charge excitations 
of the insulating system. For example the optical absorption of light with vector potential 
A ||x and finite frequency u is proportional to I]e I jx \ g)f + -^9 ^e), where \g) is 



the ground state and |e) are excited states The ground state of the insulating system 
at each cell is given by {igcr), which corresponds mainly to a hole in a 3dx2_y2 orbital (a 
Cu+^ state in agreement with spectroscopy measurements ||2^). The action of jx over this 



state gives —i{a/2)tpd{p^ a + a ) |0). The result has no overlap with states with 

«+ax/2,(T i— ax/2,(T 

one hole in Hi, but a significant overlap with states with no holes at cell i and two holes in 
a ZRS at a NN cell |(i±a yi)g2). Thus, the lowest optical excitation corresponds to energies 
near the charge-transfer gap, but slightly below it because of the attraction between the 
electron left at site i and the additional hole at a NN site. 

The dominant (resonant) contribution to the Raman intensity for incident photon energy 
uji, polarization a, and the polarization (3 of the scattered light is proportional to p5|,p^ 

J _ {f\3p\e){e\ia\g) ^2 

e Ee — Eg — Ui 



where |/) is the final state and |e) are excited intermediate states. For experimental |T9 



values of Ui and our results for the cell Hamiltonian, one can see that the excited states 
with energies Eg + uji can be those with an empty site (without holes) and another 

site with two holes: one of which is a linear combination of a 3(i^2„y2 hole and an O hig 
Wannier function, and the other can be an O hole (of symmetry a2g or 62^, see Fig. 1), 
or a non-bonding O hole. The movement of the empty site distorts the spin background, 
but, neglecting the dispersion of E^. around the local eigenenergies (eigenvalues of Hi) of 
the three two-hole states mentioned above, the intermediate spin configurations can be 
summed, and only the spin configurations of \g) and |/) enter in the expression for the 
intensity. To take into account correctly these background spin configurations, we have 
calculated the intensities from the exact wavefunctions \g) and |/) in a 4x4 cluster, in the 
Hilbert subspace of the corresponding irreducible representations. We obtain that if the 
incident frequency is near the energy required to excite an O vr hole of symmetry a2g (&2g) 
then Ia,JIb,, = 8.5 (15.7). 

The origin of the Raman excitations is similar to that described in section IV C. of Ref 
T^. For symmetry Big the O holes can also take part in the intermediate states. 

In summary, we have calculated the excitation energy and physical nature of dipole 
active E^ and novel Raman excitations of symmetries Big and A2g observed experimentally 
in insulating cuprates. The Aig, Big and Eu peaks, as well as the low-energy excitations of 
the doped systems can be well described by an effective generalized one-band model. The 
peak corresponds to a strongly bound exciton of d^y and holes. However, the amount 
of these holes in slightly doped systems is very small and does not affect the ground-state 
properties. For reasonable values of the parameters our calculated energies agree with the 
experimentally observed ones for the optical absorption peak |T^JT7[1 and Raman excitations 

We thank M. V. Klein and D. Salamon for discussions about the experimental data. Four 
of us (M.E.S., C.D.B., E.R.G., and F. L.) are supported by CONICET. A.A.A. is partially 
supported by CONICET. Partial support from Fundacion Antorchas under grant 13016/1 
is gratefully acknowledged. 



REFERENCES 

[1] CM. Varma, S. Schmitt-Rink, and E. Abrahams, Solid State Commun.621, 681 (1987); 
V.J. Emery, Phys. Rev. Lett.58, 2794 (1987). 

[2] F.C. Zhang and T.M. Rice, Phys. Rev. B 37, 3759 (1988); A.A. Ahgia, M.E. Simon and 

CD. Batista, ibzd 49, 13061 (1994). 

[3] J.J. Jefferson, H. Eskes and L.F. Feiner, Phys. Rev. B 45, 7959 (1992). 

[4] V.I. Behnicher, A.L. Chernyschev and I.V. Popovich, Phys. Rev. B 50, 13768 (1994). 

[5] M.S. Hybertsen et al, Phys. Rev. B41, 11068 (1990). 

[6] M.S. Hybertsen et al, Phys. Rev. B45, 10032 (1992). 

[7] CD. Batista and A.A. Ahgia, Phys. Rev. B 47, 8929 (1993). 

[8] H.B. Schiittler and A.J. Fedro, Phys. Rev. B 45, 7588 (1992). 

[9] M.E. Simon, M. Bahna and A.A. Ahgia, Physica C206, 297 (1993); M.E. Simon and 
A.A. Ahgia, Phys. Rev. B 48, 7471 (1993). 

[10] M.E. Simon and A.A. Ahgia, Phys. Rev. B 52, 7701 (1993). 

[11] E. Dagotto and J.Riera, Phys. Rev. Lett. 70, 682 (1993); E.S. Heeb and T.M. Rice, 
Europhys. Lett. 27, 673 (1994); CD. Batista and A.A. Ahgia, Phys. Rev. B (LT5375B). 

[12] E. Dagotto, A. Nazarenko, and M.Boninsegni, Phys. Rev. Lett. 73, 728 (1994). 

[13] CT.Chen et al, Phys. Rev. Lett. 66, 104 (1991). 

[14] B.S.Shastry and B.I.Shraiman, Phys. Rev. Lett. 65 1068 (1990). 

[15] A.V. Chubukov and D.M. Frenkel, Phys. Rev. Lett. 74 3057 (1995). 

[16] J.P.Falck et al, Phys. Rev. Lett. 69, 1109 (1992). 

[17] J.D.Perkins et al, Phys. Rev. Lett. 71, 1621 (1993). 

8 



[18] J.Lorenzana and L. Yu, Phys. Rev. Lett. 70, 861 (1993). 

[19] R.Liu et al, Phys. Rev. Lett. 71, 3709 (1993); D. Salamon et al. Phys. Rev. B51, 6617 
(1995). 

[20] M. Takigawa et al, Phys. Rev. Lett. 63, 1865 (1989); N. Niicker et al, Phys. Rev. B 
39, 6619 (1989); E. Pelegrin et al, ibid 4.7, 3354 (1993). 

[21] J.B. Grant and A.K. McMahan, Phys. Rev. Lett. 66, 488 (1991) 

[22] J.Annett et al, Phys. Rev. B 40, 2620 (1989). 

[23] L.F. Mattheis and D.R. Hamann, Phys. Rev. B 40, 2217 (1989). 

[24] J. Wagner, W. Hanke, and D.J. Scalapino, Phys. Rev. B 43, 10517 (1991). 

[25] Light Scattering in Solids II, edited by M. Cardona and G. Giintherodt ( Springer- Verlag, 
Berhn, 1982); M.V. Klein and S.B. Dierker, Phys. Rev. B 29, 4976 (1984). 

[26] G.Martmez and P.Horsch, Phys. Rev. B 44 ,317 (1991). 



9 



FIGURES 



FIGURE 1: 

Representation of the Wannier functions at a given Cu site ( denoted by a cross) constructed 
from the O 2p^ orbitals (top) and ones (bottom). 

FIGURE 2: 

Scheme of the one-hole spectral density of the extended one-band Hubbard model (right, 
full hne), the quasiparticle band for only one added hole or electron ( right, dashed line), the 
spectral density of the dxy and states (left) and the position of the different charge-transfer 
excitations (straight lines). 

FIGURE 3: 

Amplitude of the wave function of the excitons of different symmetry as a function of the 
relative coordinate between the hole and electron. The radius of each circle is proportional 
to the largest amphtude of the Hubbard one-hole creation operator of the corresponding 
cell (like \ig2) {iga\ for symmetry at any Cu site). Full and empty circles correspond to 
opposite signs. Crosses correspond to zero amplitude. For Aig we have used Vi = 0.26 to 
bind it. 
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Fig. 1 Simon et al. 




Fig. 2 Simon et ai. 



X 


O 


X 




O 


X 


• 


X • 


X 


• 


X 


o ^ 


• 


X 


o 


X • 


X 


• 


X 


• X 


X 


• 


X 


•A 


• 


X 


O 


X • 


X 


O 


e 


O ^ 


• 


X 


o 


X • 


X 


• 


X 


• X 



X 




X o • o ^ 




Fig. 3 Simon et al. 



